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Abstract 

The object of this paper is a Bayesian analysis of the 
autoregressive model X t = /?i^ t _i + £ t ,t=l,..,m and X t = 
Pi^t-i + £ t ,t=m+l,..,n where < ft, ft < 1, and e t is 
independent random variable with an exponential 
distribution with mean 8 X but later it was found that there 
was a change in the process at some point of time m which is 
reflected in the sequence after e m is changed in mean 9 2 ■ The 
issue this study focused on is at what time and which point 
the change begins to occur. The estimators of m, ft, ft and 
£^,#2 are derived from Asymmetric loss functions namely 
Linex loss & General Entropy loss functions. Both the non- 
informative and informative priors are considered. The 
effects of prior consideration on Bayes estimates of change 
point are also studied. 
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Introduction 

In applications, time series data, the common issue, 
with characteristics which fail to be compatible with 
the usual assumption of linearity and Gaussian errors, 
for a variety of reasons. Bell and Smith (1986) 
concerned with the autoregressive process AR (1) X t = 
px t _^ + £j where < /? < 1 and e t 's are i.i.d. and non- 
negative, who considered the estimation and testing 
problem for three parametric models: Gaussian, 
uniform and exponential. For large series, non 
normality may not be of importance due to the 
additive nature of filtering processes; but otherwise for 
small series. Hence, model other than Gaussian, in 
particular, the exponential, is studied. This AR (1) 
model can be used as a model for water quality 
analysis. Let the initial level of pollutant be^ ; an d 
random quantities bee 1 ,£ 2 , ■■■ Of that pollutants are 
"dumped" at regular fixed intervals into the relevant 
body of water and successive "dumping" a proportion 
(1 - P ) of the pollutant ( < /? < 1) is "washed away". 
If X t is level of pollutant at time t then ^= /J^^ + £ t 



where, £; are assumed to have an exponential 
distribution. Some of the references for Bayesian 
estimation of parameters of AR (1) as defined above 
are A. Turkmann, M. A. (1990) and M. Ibazizen, H. 
Fellage (2003). 

Apart from AR (1) with exponential errors, the 
phenomenon of change point is also observed in 
several situations in water quality analysis. It may 
happen that at some point of time instability in the 
sequence of pollutant level of a river is observed. The 
issue this study focused on is at what time and which 
point the change begins to occur, which is called 
change point inference problem. Bayesian ideas, 
playing an important role in study of such change 
point problem has been often proposed as a valid 
alternative to classical estimation procedure 

A sequence of random variables X lt ...X m ,X m+ i ...X n is 
said to have a change point at m (1 < m < n) if X t ~ 
F 1 ( x\6 x ) (i=l,2...m) and Xi ~ F 2 (x\8 2 ) (i=m+l,2...n), 
where F 1 ( x\9 1 F 2 ( x\9 2 ). The situation is in 
consideration in which and F 2 have exponential 
form, but the change point, m, is unknown. The 
problem has also been discussed within a Bayesian 
framework by Chernoff & Zacks (1964), Kander & 
Zacks (1966), A.F.M. Smith (1975), P.N. Jani & Mayuri 
Pandya (1999), Pandya, M. and Jani, P.N. (2006), 
Pandya, M. and Jadav, P. (2009) and Ebrahimi and 
Ghosh (2001). The Monograph of Broemeling and 
Tsurumi (1987) on structural changes and a survey by 
Zacks (1983) are also useful references. 

In this paper, an AR (1) model with one change point 
has been proposed, where the error distribution is 
supposed to be the changing exponential distribution. 
In section 2, a change point model related to AR (1) 
with exponential is developed. In section 3, posterior 
densities of P\,P 2 ,6\,Q 2 and m for this model are 
obtained. Bayes estimators of P\,P 2 <d l ,d 2 and m are 
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derived from asymmetric loss functions and 
symmetric loss functions in section 4. A numerical 
study to illustrate the above technique on generated 
observations is presented in section 5. In section 6, the 
sensitivity of the Bayes estimators of m when prior 
specifications deviate from the true values is studied. 
In this study, section 7 is about the simulation study 
implemented by the generation of 10,000 different 
random samples. The ending of the paper is about the 
detailed conclusion made on the study. 

f/?i*i-i + £,, i = 1,2,3 



Proposed AR (1) Model 

Let { e n , < e n , n > 1} be a random sequence having 
exponential distribution viz. 

The p. d. f . of the distribution is given by, 

/OM)= V#i.e- £iM , £ = 1.2.3 m 

= l/0 2 e~ £ i /e 2 , i = m + l, n 

Further, let {Xn, Xn, >0, n > 1 } be a sequence of random 
variables defined as, 



Xi = 



(PIXl-l 



. m 



+ £; , i = m + 1, ... .n 



(2.1) 



With X fixed constant and, < ft, ft < 1- 

In literature, there are two models for x 

Model A: x is constant, x =0 in particular; 

Model B: x has the same distribution as that of £ t . 

Here, model A for its flexibility is examined, moreover, 
the likelihood function (conditional to x ) for model B 
is exactly of the same form as that for model A. In the 

f(x Q ,x 1 x n )= (£) m . e (-^+/?i 

Since £; > 0, i = 1, 2, . . ., n. We have 



both cases, estimation will be same. 



(2.1) is the first order auto regressive process, AR (1), 
and m is unknown change point in the sequence to be 
estimated. For ft = ft =P , m = n, the model (2.1) 
reduces to the model studied by Bell and Smith (1986). 

Since, the AR (1) defined in (2.1) is Markov process, 
the joint p. d. f. of X§,X X , X m _ ,X n is given by 



(0 



. e 



<s n -s m )/8 2 



X t >ft . X t _ u and X t >0, i = 1, 2, . . ., m, X t >ft . X t _ u X t >0, 
i = m+1, . . ., n. and as a result < ft, ft < 1 . 

Then, the likelihood function giving the sample 
information is obtained 

X=(Xi,X 2 X m ,X m+ i, ...X n ) is, 

L(ft,ft, Gy, e 2 ,m\X)=^. e- A ' e \ {Q 2 Y^- m \ e~ B ^ 



9 1 ,8 2 >0, X = 



Where, 

A = S m -ftS m 

B=5 n -5 m -ft (S n - S m ) 

Bayes Estimation 



(2.2) 



The ML methods as well as other classical approaches 
are based only on the empirical information provided 
by the data. However, when there is some technical 
knowledge on the parameters of the distribution 
available, a Bayes procedure seems to an attractive 
inferential method. The Bayes procedure is based on a 
posterior density, say, 5 (ft, ft, 6\, 02., m \ X), which is 
proportional to the product of the likelihood function 



.02. 

Sic - 2~]?=i 

^(ft>ft> &i> 2 ,m | X), with a prior joint density, say, 
5 (ft- ft- &i> &2> m I K) representing uncertainty on the 
parameters values. 

9(Pi,p2,0 1 ,9 2 ,m\X) = 

LQhjhA&W | X) .g{Pl,Pl,<h,th,m I jQ 

YSa=lSp j/p 2 ! lh ! 02 L{Pl.P2A,Oz,m \X).g{p 1 ,p 2 ,<h,th,m \X)dp 1 dp 2 de 1 d% 

Using Informative Priors for pi and ft, Oi and 9i 

The first step in a Bayes analysis is the choice of the 
prior density on the distribution parameters. When 
technical information about the mechanism of process 
under consideration is accessible, then this information 
should be converted into a degree of belief on the 
distribution parameters (or function thereof). This 
problem is greatly aggravated when some 'physical' 
meaning can be attached to these parameters. 

It is also supposed that some information on these 
parameters is available, and that this technical 
knowledge can be given in terms of prior mean values 
(ii, |i2. Independent beta priors on ft and j}2 with 
respective means jul, ju2 and common standard 
deviation a viz are supposed. 

5(ft) = [ft ai " 1 (l- ft)" 1 " 1 ] I P(ai,&i), 
<ft<l, a x >0, b x >0 
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5(i5 2 ) = [/?2 2 - 1 (l- ft)* 2 " 1 ] I (3(a 2 ,& 2 ), 
< /? 2 < 1, a 2 > 0, b 2 > 



andi?!, 6*2- 



If the prior information is given in terms of prior 



g(m) = 



n — 1 



means \xi and a, then the hyper parameters can be For simplicity, it is assumed that, a priori, 

ft, ft, ft, #2 and m are independently distributed. The 
joint prior density of ft, ft, ft, ft and m say, 



0i(ft,ft,ft, ft,m) is, 

»<A.ft «. ft.™) =t , g^z£i2^|p!z£^ 

Where, /q = 



Pla 1 b 1 )P (a 2 6 2 )(n-l) 



obtained by solving 

Of = o-j -1 [C 1 ~ Aii)i"i 2 ~ i = 1/2- 
&i = juf 1 (1 - Hi)ai i = 1,2. (3.1) 
Let the joint prior density of ft and 9 2 be 

As in Broemeling et al.(1987), the marginal prior Joint posterior density of ft, /? 2 , ft, ft and m say, 
distribution of m is supposed to be discrete uniform fl'ivPi.ft. ft. ft, m l *) is, 

over the set { 1, 2, n - 1} and ndependent of f>l, j}2 

9i (ft . ft . ft. ft. ™ I X) = h 0[ ( m+1 W* ^T^" 1 e -B/ft 

ft^-Hi- ft)' 1 " 1 p^ 2 -Ki- P 2 ) b2 - 1 /h(x) , 

hi(x) is the marginal density of X_ given as, where, 
h(x) = k 1 Z n m= \T 1 (m) , (3.2) 

Ti(m) = = Tm r(n- m). f(-)(5 m )" m /IppeZ a^m, - ^ , 1+ a v ^ ,1]} 

((f) [(5„ - 5 m )]-("" m ) J 4 ppe Z 1 F 1 [a 2 „n-m,-fe 2 ,l + a 2 ,^|^,l]] , (3.3) 



And, 
Note-1 



Appel 1 F 1 (a,b 1 ,b 2 ,c,x,y) is Appel Hypergeometric 
function defined as 



Appel 1 F 1 (a,b 1 ,b 2 ,c,x,y) = — ^ f u^Hl - u) c_a_1 (l - uxY^O- - uyY^du, 

ia i\c — a) J 



r<a r[c - a) 
Real[a] > 0, Real (c-a) > 0. 

Integrating g x (ft.ft, ft, ft,m | *) with (ft, ft) and 
(ft, ft), leads to the posterior distribution of change 
point m. 



T^m) is as given in (3.3). 



(3.4) 



Marginal posterior densities of ft and ft say, ^(ft I X) 
and g 2 (ft IX) are as, 



Marginal posterior density of change point m is given 

g 1 (p 1 I X) = r m -=\ Jo 1 J" " Jo" 9i (Pi.fo. ft, (h,m I X) .de 1 de 2 dp 2 
ft^-Hi - ft)" 1 " 1 rjl£ ^P 1 {{^} Ks n - s m )]~^ 

.Appel [ a 2 „n - m, - 62, 1 + ^.^Tr- 1]] / hiQO 

5l (ft 1 jo -*i z sri ft 02 " 1 (i - ft)" 2 " 1 ^ !! ^ 2 {(£) (5 m )- m 

Appel 1 F 1 [a 1 ,m,-b 1 , 1+ ^,^,1]]/ ^Qf) 
Marginal posterior densities of and 2f say ^ IX) and gi(9 2 \ X_) are given as, 

9i{o 1 \x) = r m -=\ Jo 1 Jo 1 j7 <?! (ft, ft, ft, ft.™ 1 r ) de 2 dp l dp 2 



(3.5) 



(3.6) 
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^lETA 0-, im+l) ^-Jf 1 " d Pl nn- m){ (f) [(S„ - S m )]-(— ) 



a 2 ; 



A VV el x F x \a 2 ,,n-m-b 2 ,\ + a 2 ;f^,\\l hQQ (3-7) 



And 



5i (<% i x ) = sr=\ Jo 1 Jo 1 foai (Pi.fo, Ox, e 2 ,m i J0de 1 dp 1 dB 2 

= ^a-i ^ ( -v m / /° 2 "y d fe {(i) ( 5 m 

Appel 1 F 1 [a i ,m,-b 1 , 1 + ,1] ]/ hiQQ (3.8) 

Using Non Informative Priors for pi, $z, 8i and 82 variables, then the non-informative prior is, 

The non-informative prior is a density which adds no 9i (Pi, Pi, #l. 9 2 ,m ) K 9 (n _ i:) 

information to that contained in the empirical data. If 

., , , . c n n n n Joint posterior density of 81, 82, Bi, 82 and m say, 

there is no available information on ft, ft, 0\, &i, m r : 11 3 

which are assumed to be a priori independent random 92 02> &i> &2> m \ K) ^ s given by 

= J_g-(m+l) e s m+ p 1 st n /9 l ^-(n-m+1) e -S n +S m +/? 2 (S;-S^)/<% /^(X) 

Where, ft 2 (x) is the marginal density of x given by, 

h 2 {x)=^ ] Y l n m -} 1 T 2 (m) (3.9) 

Where, 

T 2 (m) = rmr(n - m) \- ^ + ( f m ~ sSl ^ f- (5 " ~ 5 ™ >~ ( " ~ m )+ ^ zga -Jn+s^ n ~ m ^ (3 10) 

V ; V V L (Sm)(m) V (n-m)(S*-S^) /J v ; 

Marginal posterior density of change point m is given by 

n-l ~ =0 1 1 

52 (m I X) = Yjjjjj 92 (ft,/?2 ' ^' 6 ' 2 ' m ' ^ de i d0 2dP 1 dp 2 

m=l 

= Tj^/SrArzCm) (3.11) 
Marginal posterior density of ^and 9 2 are given as, 

9 2((h 1 x) =j- a-i *r (m+1) ■ ( - 1) 

f (*n - S m )-(" ~"0 + (5„ -S m -S* )-(" - m ) \ r x 

I (n-m)(s;-5«) J/M*J (3-12) 

9l(0 2 1 X) 

(s„- s, 



£ SSI**** ^ I-^IF]/ ^) (3.13) 



Marginal posterior density of ft and ft are given as, 



52 (Pi I £ Sr=\ rmr(n- m) 



IV (n-m)(5*-5*) JJ/M*J (3-14) 



(n-m)(5*-5^) 

o« 1 S ■ ^ K='. r„r(„ - ^ __l_^ -'~> 

k- 5 ""';!",; 5 ;'"" )] Mx) (3.i5) 

Remark 1: For n=m, ft = ft, Q x = 9 2 the equations (3.5), densities of /? and 9 of AR(1) process without change 
(3.7) and (3.12), (3.14) reduce to the marginal posterior point. 
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Bayes Estimates under Symmetric and 
Asymmetric Loss Functions 

The Bayes estimate of a generic parameter (or function 
thereof) a based on a Squared Error Loss (SEL) 
function Li (a, d)= (a-d) 2 , where d is decision rule to 
estimate a which is the posterior mean. As a 
consequence, the SEL function is relative to an integer 
parameter, 



Li'(m, v) °° (m-v) 



m, v=0,l,2,. 



Hence, the Bayesian estimate of an integer-valued 
parameter under the SEL function Li'(m, v) is no 
longer the posterior mean and can be obtained by 
numerically minimizing the corresponding posterior 
loss. Generally, such a Bayesian estimate is equal to the 
nearest integer value to the posterior mean. So, the 
nearest value to the posterior mean is regarded as 
Bayes Estimate. The Bayes estimators of m under SEL 
are the nearest integer values to (4.1) and (4.2), 

m* = JZrAmTdm) / ET^C"*) 



m = 

Zl- = \mT 2 (m) /Zl- = \T 2 (m) 



(4.1) 



(4.2) 



Other Bayes estimators of a based on the loss functions 
Lz (a, d) = I a - d I 

L 3 (a, d) = 



if | a - d | < e, e>Q 
otherwise 



is the posterior median and posterior mode, 
respectively. 

Asymmetric Loss Function 

The Loss function L (a, d) provides a measure of the 
financial consequences arising from a wrong decision 



rule d to estimate an unknown quantity a. The choice 
of the appropriate loss function is just dependent on 
financial considerations, but independent upon the 
used estimation procedure. In this section, Bayes 
estimator of change point m is derived from different 
asymmetric loss function using both prior 
considerations explained in section 3.1 and 3.2. A 
useful asymmetric loss, known as the Linex loss 
function was introduced by Varian (1975). Under the 
assumption that the minimal loss at d, the Linex loss 
function can be expressed as, 

L4 (a, d) = exp. [qi (d- a)] - qi (d - a) - 1, qi * 0. 

The sign of the shape parameter qi reflects the 
deviation of the asymmetry, qi > if over estimation is 
more serious than under estimation, and vice-versa, 
and the magnitude of qi reflects the degree of 
asymmetry. 

Minimizing expected loss function Em [L4 (m, d)] and 
using posterior distribution (3.4) and (3.11), the Bayes 
estimate of m using Linex loss function is obtained by 
means of the nearest integer value to (4.3) under 
informative and non-informative prior , say, say m* L 
m* L * . 

ml=- %! /»ET=i e-^riCmVXTiiTUm)], 

mV = ~ V<h 'n ESTii e ~ aimT 2 ("0/ETii T 2 (m)], (4.3) 

Where Ti (m) and 7 2 (m)are as given in (3.3) and 
(3.10). 

Minimizing expected loss function Em [L4 (pi, d)] and 
using posterior distributions (3.5) as well as (3.6), the 
Bayes estimators of jSi and jS2 are obtained using Linex 
loss function as, 



p{ L = In [£e"«''] 



= -i- In [XriJo 1 ' 1 " 1 " 1 ^! 1 " 16 "" 1 ' 1 dp, rmr{n- m){( ^) [{S n - S m )]-<»-> 



a 2 

Appel 1 F 1 [a 2 ,,n-m,-b 2 ,l + a 2 ,f^,l]}/h 1 (,)Q 



(4.4) 



n-l 

P* 2L ==-^ln [V f 



V2 2 - 1 (l-p 2 ) b2 - 1 -e 



-11P2 



B 



n—m + 1 



dp 2 rmr(n - m) { (S m )" 



Appel 1 F 1 [a 1 ,m,-b 1 , 1+ ,1]} / h^X) 

Where Appel ^[a^m,- b x , l+a,,^-,!]} and 

1] are same 



(4.5) 



Appel x ¥ x \ a 2 ,,n-m,-b 2 ,l + a 2 ,j-^r 
as in note 1 and h x (X) is as in (3.2). 



Another loss function, called General Entropy loss 
function (GEL), proposed by Calabria and Pulcini 
(1996) is given by, 



57 



www. srl-j ournal . org 



Statistics Research Letters (SRL) Volume 2 Issue 2, May 2013 



Ls(a, d) = (d/ci) qi - q 3 l n ( d / a) - 1, using General Entropy loss function by means of the 

.... nearest integer value to (4.6) and (4.7) under 

Whose minimum occurs at d = a, minimizing ° . . . 

_ rT . ,., , . , informative and non-mformative prior, say m P and 

expectation b [Ls (m, d)J and using posterior density „ 

gi(m | x ), i = 1,2. The Bayes estimate of m is achieved 

m\ = [E, [m-^]]~ 1/q3 = Er=\ m-^Tl (m)/^ 1 ! 7i (m)Y l/q \ (4.6) 

m? = [m-^]]~ 1/q3 = E^ii m-«r 2 (m)/^ 1 ! ^2 (m)]" 1 ^, (4.7) 

Putting q3 = - 1 in (4.6) and (4.7), Bayes estimate, using posterior distributions (3.5) and (3.6), we get the 
posterior mean of m is acquired. Bayes estimates of pi using General Entropy loss 

Minimizing expected loss function E [Ls (pi, d)] and function respectively as, 

P* = METiirm r(n - m)-±- (s m r m ... * f [(s n - s m )]~^ 



Appel 



a 1 -q 3 m,-b 1 , 1 + a t - q 3 ,,-^-, 1 



J 4ppeZ 1 F 1 [a 2 „n-m,-b 2 , 1 + 02,^^,1]] Mi®} V<?3 (4.8) 
^=feiSr=\rm r(n - m) [(5 n - S m )]-("" m ) f{(S m )" m 

^* ^* 

/IppeZ iFi[a 2 -q 3 ,,n-m, - b 2 ,\ + a 2 - q 3 , " _ m , 1 ]} 

Appel 1 F 1 [a 1 ,m,-b 1 , 1 + a v ^- /h^X)}' 1 ^ 3 (4.9) 

Where ,4ppeZ t F t [ m, - b t , 1+ a u -f- ,1]] and using posterior distributions (3.14) and (3.15), we 

s*-s* m get the Bayes estimate (3£ of pi, i = 1, 2 using non 

Appel iFiT a 2 ,n - m, — b 2 ,,l + a 2> — — -,11 are same K . c .. . , ^ , c , T 

^ i il /' z- Jn _ Jm - j informative priors, under General Entropy Loss 

as in note 1 and h x (X) is as in (3.2). function is given as, 

Minimizing expected loss function E [Ls (pi, d)] and 

Pie - l^I^m=i im K (n _ m)(5 ._^) )\ 

i/q 3 



(S m r m 2 F 1 [l-q 3 ,m,2- q 3 .^](l-?3)- 1 /h2fe)} (4-10) 

B „ f 1 Vrr ^ S- m + (S m - S^)- m 
2E (n-lL C%)(m) 



(5 B - 5 m )-("- +1 ) jPitl-,3. n-m, 2- q3 ,^ ](1 _ q3 )-i/ h 2 (x) (4.11) 

Where ft 2 (^) is same as in (3.9). using posterior distributions (3.7) as well as (3. 12), we 

get the Bayes estimate 0[ E and ff{ E of 0i using General 
Minimizing expected loss function E [Ls (9i. d)l and c , T c .. 

b r Entropy Loss function as, 



<?* = fci y r(m + q 3 ) T(n - m) 1 i(5 m )-( m+ ^)[(5 n - 5 m )]-("" m ) 



m=l 



ilppeZ 1 F 1 fa 1 ,m,-6 1 , l+cn.^-.ll] ^peZ ^ a 2 ,n - m, - fe 2 „ 1 + a 2 , S f^, l]/ ^(x)}"" (4.12) 

- Z m=1 T(n - m) - [ (n- m)( s«- S «) J 

V - Sm > 0,s m > 0,m + q 3 >0, 
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Minimizing expected loss function E [Ls (02, d)] and the Bayes estimate 9* 2E and 9* 2E of 02 using General 
using posterior distributions (3.8) and (3. 13), we get Entropy Loss function as, 

&2E = {fci HTii rmTin - m + q 3 )(f ) [tf„ - S m )Y^^ (±) (S m )"( m+1 > 

c* _ c* 

■Jtj ^ hi 



Appel t F t 



/?* — 



a 2 ,n-m + q 3 ,- b 2 ,,l+ a 2l - _ 
Appel jFi [ ai , m, - b x , 1 + a x ,^ ,1]} /h x (x)}~ 
r(-l + n - m + q 3 ) 



1/93 



(4.14) 



n-lL 



0>n ) 



1+n— m+173 



+ 



( 1 - - ) 



-1+n—m +^3" 



c-m , (c* _ c* \- 



Mi)} 



-1/93 



(S^)(m) 

> 0, s m > 0, m + q 3 >0, 



(4.15) 



Remark 2: For n=m, /?! = /? 2 , 6*i = #2 the equations (4.4), 
(4.8), (4.10) and (4.12), (4.13) reduce the Bayes 
estimates under asymmetric loss function of /? and 8 of 
AR(1) process without change point under asymmetric 
loss functions. 

Numerical Study 

Let us consider AR (1) model as 



X l = { 



i = 1,2, ...12 
i = 13,14,... 20. 



Where, £ ; 's are independently distributed exponential 
distributions given in (2.1) with 

8i=0.6, 82=1. 20 random observations have been 
generated from proposed AR (1) model given in (2.2). 
The first twelve observations are from exponential 
distribution with 8i=0.6 and next eight are from 
exponential distribution with 82=1. pi and p2 
themselves are random observations from beta 



distributions with prior means |ii = 0.1, 

TABLE 1 GENERATED OBSERVATIONS FROM PROPOSED AR (1) MODEL. 



U2 = 0.8 and common standard deviation a = 0.1 
respectively, resulting in ai = 0.001, 

b =0.09, a2=0.48 and b2= 0.12. These observations are 
given in Table-1. 

Posterior mean of m, 01, 91, jil, and jil and the 
posterior median of m have been calculated.. Posterior 
mode appears to be a bad estimator of m. For a 
comparative purpose point of view, estimators under 
the non informative prior are also calculated. The 
results are shown in Table-2. The Bayes estimates 
m* L ,m* E of m, 6\ E , 6* 2E of 9\ and 6 2 , P[ E , P 2E of jil and 
fa respectively for the data given in Table-1 have been 
computed. As well as the Bayes estimates under the 
non informative prior and asymmetric loss functions 
which are also calculated, from which the result shown 
in Table-3 reveal that m\, m*E, fi{ E , fS 2E , 9\ E and 6* 2E 
are robust with respect to the change in the shape 
parameter of GE loss function. 
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1.21 


0.44 


I 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


6, 


0.90 


0.99 


1.07 


0.15 


2.05 


0.04 


0.14 


0.86 


1.12 


0.32 


I 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


Xi 


0.38 


0.41 


0.86 


1.68 


3.01 


2.53 


4.36 


5.53 


4.57 


4.61 


I 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


e. 


0.33 


0.43 


0.48 


0.99 


1.67 


0.12 


2.33 


2.04 


0.15 


0.95 
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TABLE 2 THE VALUES OF BAYES ESTIMATES OF CHANGE POINT 



Prior Density 


Bayes estimates of change point 


Bayes estimates of 
auto correlation 
coefficient. 


Bayes estimates of 
Exponential parameters 


Posterior Median 


Posterior Mean 


P» 


6 2 


9i 


02 


Informative 


12.02 


12 


0.1 


0.8 


0.6 


0.95 


Non informative 


12.22 


11 


0.12 


0.82 


0.63 


1.03 



TABLE 3 THE BAYES ESTIMATES USING ASYMMETRIC LOSS FUNCTIONS. 



Prior Density 


Shape 


Bayes estimates of change 


Bayes estimates under General Entropy 




parameter 


point 




Loss 






qi 


q 3 


m' L 


m* E 


P*2E 


A* £ 






Informative 


0.8 


0.8 


12 


12 


0.83 


0.13 


0.13 


0.95 




1.2 


1.2 


12 


12 


0.82 


0.12 


0.12 


0.94 




1.5 


1.5 


12 


12 


0.81 


0.11 


0.11 


0.91 


Non 






m L 


m E 


m 


m 




$IE 


informative 


0.8 


0.8 


13 


13 


0.86 


0.16 


0.68 


1.4 




1.2 


1.2 


13 


13 


0.85 


0.15 


0.66 


1.3 




1.5 


1.5 


13 


13 


0.83 


0.13 


0.65 


1.2 



TABLE-4 POSTERIOR MEAN M* FOR THE DATA GIVEN IN TABLE-1. 





U2 


m* 


m*E 


0.1 


0.6 


12 


12 


0.07 


0.8 


12 


12 


0.2 


0.4 


12 


12 



Sensitivity of Bayes Estimates 

In this section, the studied sensitivity of the Bayes 
estimates is obtained from section 4 with respect to 

change in the prior of the parameter. The means jU^ , 

jj,2 an d standard deviation o of beta prior have been 
used as prior information in computing the parameters 
ai, bi, ai, hi of the prior. Following Calabria and Pulcini 
(1996), it is also assumed that the prior information 
should be correct if the true value of Pi (P2) is close to 
prior mean lu(lu) and is assumed to be wrong if Pi (P2 ) 
is far from ni(ua). We have computed m* and m* E using 
(4.1) and (4.6) for the data given in Table-1 with 
common value of =0.1. For q3 =0.9, considering 

different values of ( jJ,y , /J, 2 ) an ^ result are shown in 
Table-4. 

The results shown in Table-4 lead to conclusion that m* 
and m*E are robust with respect to the correct choice of 
the prior density of Pi (P2) and a wrong choice of the 
prior density of Pi (P2). Moreover, they are also robust 
with respect to the change in the shape parameter of 
GE loss function. 

Simulation Study 

In section 4, we have obtained Bayes estimates of m on 



the basis of the generated data given in Table-1 for 
given values of parameters. To justify the results, we 
have generated 10,000 different random samples with 
m=12, n=20, 6i=1.0, 02=2.0, |3i= 0.4, (3 2 = 0.6 and 
obtained the frequency distributions of posterior mean, 
median of m, mL*, mE* with the correct prior 
consideration. The results are shown in Table 5. We 
also obtain the frequency distributions of Bayes 
estimates of autoregressive coefficients given in section 
4 with the both prior considerations and the results are 
shown in Table 6 and table 7. The value of shape 
parameter of the general entropy loss and Linex loss is 
taken as 0.1. 

We also simulate several samples from AR (1) model 
explained section 2 with m=12, n=20, 61 =1.0, 0.6, 0.7; 62 
=2.0, 0.8, 0.9 and (3i=0.2, 0.4, 0.5; (3 2 =0.4, 0.6, 0.7. For 
each 0i, 02, pi and p2, and Bayes estimators of change 
point m and autoregressive coefficients pi and p2 using 
q3= 0.9 has been computed for different prior means \ii 
and L12, Which lead to the same conclusion as for single 
sample, that the Bayes estimators posterior mean of m, 
and m*E are robust with respect to the correct choice of 
the prior specifications on Pi (P2) and wrong choice of 
the prior specifications on p2 (Pi) 
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Table 5: Frequency distributions of the Bayes estimates of the change point 



Bayes estimate 


% Frequency for 




01-10 


11-13 


14-20 


Posterior mean 


18 


70 


12 


Posterior median 


10 


74 


16 


Posterior mode 


12 


73 


15 


mL* 


13 


77 


10 


mE* 


14 


78 


8 



Table 6: FREQUENCY distributions of the bayes estimates of autoregressive coefficients Bl and b2 using general entropy loss function 

(INFORMATIVE) 



Bayes estimate 


% Frequency for 




0.1-0.3 


0.3-0.5 


0.5-0.7 


0.7-0.9 




13 


82 


02 


03 


P1e 


11 


01 


84 


04 


Table 7: frequency distributions of the bayes estimates of (non-informative) 


Bayes estimate 


% Frequency for 




0.1-0.3 


0.3-0.5 


0.5-0.7 


0.7-0.9 




13 


78 


06 


03 




11 


01 


74 


14 



Conclusions 

Our numerical study shows that the Bayes estimators 
posterior mean of m, and m*E are robust with respect to 
the different choice of the prior specifications on Pi, P2 
and are sensitive in case prior specifications on both (3i, 
P2 deviate simultaneously from the true values. 
Numerical study also shows that posterior mean of m 
is sensitive when prior specifications on both (3i, P2 
deviate simultaneously from the true values. 
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